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ABSTRACT 

Formation of sub-galactic halos is suppressed in warm dark matter (WDM) model due to thermal 
motion of WDM particles. This may provide a natural resolution to some puzzles in standard cold dark 
matter (CDM) theory such as the cusped density profiles of virialized dark halos and the overabundance 
of low mass satellites. One of the observational tests of the WDM model is to measure the gravitationally 
lensed images of distant quasars below sub-arcsecond scales. In this Letter, we report a comparison of 
the lensing probabilities of multiple images between CDM and WDM models using a singular isothermal 
sphere model for the mass density profiles of dark halos and the Press-Schechter mass function for their 
distribution and cosmic evolution. It is shown that the differential probability of multiple images with 
small angular separations down to ~ 10 milliarcseconds should allow one to set useful constraints on the 
WDM particle mass. We discuss briefly the feasibility and uncertainties of this method in future radio 
surveys (e.g. VLBI) for gravitational lensing. 

Subject headings: cosmology: theory — dark matter - galaxies: formation — gravitational lensing 



1. INTRODUCTION 

While cosmological models with a mixture of roughly 
35% cold dark matter (CDM) and 65% vacuum energy 
have proved remarkably successful at explaining the ori- 
gin and evolution of cosmic structures on large-scales (> 1 
Mpc), there are a number of independent observations on 
small scales that seem to be in conflict with the predic- 
tions of standard CDM models such as the cusped density 
profiles of virialized dark halos and the overabundance of 
satellite structures. Among many other mechanisms sug- 
gested for resolving these discrepancies, the simplest ap- 
proach is to smear out small scale power by free-steaming 
due to thermal motion of particles. Warm dark matter 
(WDM) is thus invoked and has recently received a lot of 
attention in literature (e.g. White & Croft 2000; Colin, 
Avila- Reese & Valenzuela 2000; Sommer-Larsen & Dolgov 
2000; Bode, Ostriker & Turok 2000b; references therein). 

In the conventional scenario of structure formation by 
WDM, perturbations on small mass scales Mf « 4 x 
lO n M /i 2 fi x (.Rf/O.l Mpc) 3 are damped, relative to CDM 
models, where fix is the WDM density parameter and 
Ri denotes the characteristic free-streaming length. As 
a result, it suppresses the formation of small halos with 
M < Mf at high rcdshift and meanwhile reduces the over- 
all number of low mass halos in the universe, depending 
sensitively on the WDM particle mass mx- Observational 
tests of this scenario and hence robust constraints on mx 
are restricted to the studies of less massive systems with 
M < Mf. The existing limits in literature are based on 
the Gunn-Peterson effect and the Ly-a forest, which set 
mx > 750 eV (Narayanan et al. 2000). In this Letter, 
we would like to examine whether the gravitational lens- 
ing of high-redshift quasars by foreground dark halos can 
be used to distinguish CDM and WDM models or place 
any meaningful constraints on mx- This arises because 
the probability of gravitational lensing would be lower in 
WDM model, as a result of the suppression of formation of 



small mass halos on galactic and sub-galactic scales, than 
in CDM model. For a massive galaxy like the Milky Way 
as a lens at cosmological distance, the Einstein radius of 
a background quasar is typically 1 arcsecond. So, it is 
expected that the difference in lensing properties between 
WDM and CDM models would occur below sub-arcsecond 
scales. We would like to demonstrate quantitatively to 
what angular separations the two models become to be 
distinguishable and how sensitively their lensing probabil- 
ities depend on the WDM particle mass. 



2. LENSING PROBABILITY 

For simplicity, we use a singular isothermal sphere to 
model the mass distribution inside a virialized halo so 
that the total mass M within radius r is simply M = 
2a 2 r/G. We truncate the mass profile at the virial ra- 
dius r v ; r defined by M v ; r = 47rA c( o c r 3 , ir /3, where p c = 
3H 2 (z)/8itG is the critical mass density of the universe, 
H(z) = 100 h E(z) km s" 1 Mpc" 1 is the Hubble constant, 
E(z) = f2 M (l + z) 3 + n K (l + z) 2 + n A , n K denotes the 
curvature term, and A c represents the overdensity of dark 
matter with respect to the critical density p c at the col- 
lapse redshift z. For a dark halo at z as a lens, the multiple 
images of a distant quasar at z s with an angular separation 
of A8 — 26>e will occur if it happens to lie in the Einstein 
radius of the dark halo: (9e = ■in(a v /c) 2 (D ( i LS /D s ), in which 
Dds and D s are the angular diameter distances from the 
lens and observer to the source, respectively. The lensing 
cross section for the multiple images is simply 7r# E . The 
probability that a quasar at z s will have multiple images 
with an angular separation greater than AO due to all the 
foreground dark halos is thus 
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where dV is the comoving volume element, and M min is 
related to AO through 

M min = 1.50 x 1O 12 M hr 1 




We use the Press-Schechter (PS) mass function for the 
abundance and evolution of virialized dark halos: 

dN flj^5_c(z)do_ ( 5 c (z) 2 \ 

dMdV \irM ^ dM eXP { 2a 2 )' [ ) 

where p is the present-day average mass density of the 
universe, S c (z) is the linear over-density of perturbations 
that collapsed and virialized at redshift z, and er is the 
variance of the mass density fluctuation in sphere of mass 
M = (47r/3)pi? 3 . It is convenient to express 6 c (z) as 
S c (z) = Sg (1+z) [g(0)/<?(z)], where g(z) is the linear growth 
factor, for which we take the approximate formula by Car- 
roll, Press & Turner (1992), and the quantity So has a weak 
dependence on Qm- S n = 1.6866[1 + 0.01256 log^M^)] for 
a flat universe with nonzero cosmological constant (Math- 
iesen & Evrard 1998). Once a power spectrum, P(k), is 
specified, the mass variance becomes 

1 f°° 

a 2 {M) = — / k 2 P(k)W 2 (kR)dk, (4) 

in which W{x) = 3(sinx — xcosx)/x 3 is the Fourier repre- 
sentation of the window function. We adopt the following 
parameterization of WDM power spectrum 

P(k) = Ak n T 2 (k)T 2 BM (k), (5) 

where n is the primordial power spectrum and is assumed 
to be the Harrison- Zeldovich case n = 1, TcDu(k) is 
the transfer function of adiabatic CDM model for which 
we use the fit given by Bardeen et al. (1986), and 
T x (fc) acts as the CDM to WDM " transfer function" 
for which we adopt the form provided by Bode et al. 
(2000b): T x (fc) = [1 + {ak) 2v ]~ 5 / v ', where v = 1.2 and 
a is related to the WDM particle mass mx through 
a = 0.048(fi x /0.4) 015 (V0.65) 1 - 3 (keV/m x ) 115 . The am- 
plitude A is determined by equation (4) using the rms mass 
fluctuation on an 8 ft -1 Mpc scale, a^. 

We work with a flat cosmological model where Om = 
O x + n h = 0.3, n h = 0.03, fl A = 0.7, and h = 0.68, 
which fixes the shape parameter T — Q,Mhexp[— f2t>(l + 
\/2/i/Om)] = 0.176. The normalization parameter is taken 
from the calibration of cluster abundance, er 8 = 0.85. We 
compute the lensing probabilities for a quasar at redshift 
z s = 1, 2 and 3, respectively, and for two choices of the 
WDM particle mass mx — 350 eV and 1 keV. As a compar- 
ison, we also give the corresponding lensing probabilities 
for a CDM model with Oqdm = 0.27. In order to highlight 
their differences, we display in Figure 1 the normalized dif- 
ferential probabilities [l/P(A6)][dP(A6)/dA6] instead of 
the total probabilities P(A9). For the latter, the typical 
value for a quasar at z s = 2 with AO > 1" is nearly the 
same for all the three models (two WDM models and the 
standard CDM model): P(> 1") w 1 x 10~ 3 . This re- 
sult is consistent with the pioneering statistical study of 



gravitational lensing by Turner, Ostriker & Gott (1984). 
Noticeable departure of the lensing probability in WDM 
models from that in CDM model occurs only at small an- 
gular separations AO < 1" and AO < 0."1 for m x > 350 
eV and mx > 1 keV, respectively. It is likely that one 
needs to reach an even smaller angular separation of ~ 10 
milliarcseconds in order to make an effective distinction 
between the WDM and CDM models for a realistic WDM 
particle mass of mx ~ 1 keV. We have applied the lens- 
ing probability P{A0) to the quasar luminosity function 
determined from the 2dF QSO survey over absolute mag- 
nitudes -26 < M B < -23 and redshifts 0.35 < z s < 2.3 
(Boyle et al. 2000). Here we would rather take a less 
vigorous approach to the problem by neglecting the mag- 
nification effect of the double images. It is shown that at 
a limiting magnitude of me = 20, the fractions of small 
angular separation events with 0."001 < AO < 0."01 in 
the total number of doubly-imaged quasars are approxi- 
mately 1/100, 4/1000 and 1/1000 for the CDM model and 
the WDM models with mx = 1 keV and mx = 350 eV, 
respectively. Indeed, observational tests of these fractions 
on milliarcsecond scales turn out to be difficult and chal- 
lenging. 

3. DISCUSSION AND CONCLUSIONS 

In the CDM scenario of structure formation, many small 
galaxies would form, relative to what are predicted by the 
WDM model. Consequently, the standard CDM model 
would predict more multiple, small angular separation im- 
ages of the gravitationally lensed distant quasars than 
WDM model does. This unique property can be used as 
an effective tool to distinguish these two prevailing mod- 
els. It has been shown that the significant difference in 
the distributions of image angular separations between 
WDM and CDM models occurs at small angular sepa- 
rations AO <~ 0.01" for the WDM particle mass mx ~ 1 
keV. Although there are no data available at present to test 
our prediction, these angle ranges should be accessible by 
current radio observations such as VLBI. Recall that sys- 
tematic radio surveys for strong gravitational lensing have 
proved to be very successful at finding small separation 
images in the range 0."3-0."6 (e.g. Augusto, Wilkinson & 
Browne 1998; Augusto & Wilkinson 2001; and references 
therein). It is hoped that future multiple image surveys 
down to ~ 10 milliarcseconds would allow one to place 
useful constraints on the WDM particle mass. 

We now briefly discuss several major uncertainties in 
the present predictions: Firstly, we have used a simple 
singular isothermal sphere model for the mass distribu- 
tions of dark halos. An alternative is the universal density 
profile suggested by high-resolution A-body simulations 
(Navarro, Frenk & White 1995; NFW), which should be 
valid for the virialized dark halos as small as 10 9 M Q . Li & 
Ostriker (2000) have recently compared the lensing prob- 
abilities P(A0) produced by a singular isothermal sphere 
model and the NFW profile using roughly the same ap- 
proach as ours. It turns out that the NFW profile results 
in a significantly small value of P(A0) ~ 10~ 6 , insensitive 
to the image separation AO. We have recalculated the lens- 
ing probabilities under CDM and WDM models following 
the approximate treatment of Li & Ostriker (2000) for the 
image separation by NFW profile, AO w 20q where 0o is 
the solution to the lensing equation with the zero align- 
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ment parameter. Unlike Li & Ostriker (2000) who took a 
constant halo concentration c for massive halos from nu- 
merical simulations, we fix the characteristic density S c 
and halo concentration c for each halo M using the S c -M 
relation via the collapse redshift z co \\ proposed by Navarro, 
Frenk & White (1997). Nevertheless, we have reached a 
conclusion essentially similar to Li & Ostriker (2000). For 
example, the lensing probabilities for a distant quasar at 
z s = 3 with AO > 0.001" - 1" are P(A9) w 3.7 x 10" 6 , 
3.6 x 10~ 6 and 2.9 x 10~ 6 for the CDM model and the 
WDM models with mx = 1 keV and mx = 350 eV, respec- 
tively. The lack of small separation events below arcsecond 
scales will invalidate our original intention that the study 
of doubly-imaged quasars may help to distinguish between 
CDM and WDM models. However, we note that the lens- 
ing probability produced by NFW profile is unreasonably 
small (~ 10~ 6 ). This implies that the NFW profile may be 
inappropriate to the description of the mass distributions 
in the inner regions of dark halos where the multiple im- 
ages would occur. Indeed, the central mass density of an 
NFW profile is too shallow to act as an effective lens. The 
absence of detectable odd images in the known lens sys- 
tems also lends support to the existence of a steeper inner 
density profile in galaxies, p cx r -2 (Rusin & Ma 2001). 
In a word, the singular isothermal sphere model may be 
a reasonable approximation for the inner mass profiles of 
galactic and sub-galactic halos. 

Secondly, the PS mass function may become inaccurate 
for small halos. High resolution TV-Body simulations have 
shown that the PS approximation predicts too many small 
halos at low redshifts (e.g. Bode et al. 2000a; Jenkins et 
al. 2001). This will lead to an overestimate of the total 
lensing probability for multiple images, P(A6), although 
the differential lensing probability normalized by P(A9), 
[l/P(A6)][dP(A6)/dA9], is probably less affected. In or- 
der to demonstrate the uncertainty, we have used the mod- 
ified PS mass function by Sheth & Tormen (1999) and 
the best fit parameters by Jenkins et al. (2001). It is 
found that the total lensing probabilities for small separa- 
tion events in both CDM and WDM models are slightly 



reduced as compared with those predicted by the standard 
PS formalism. For example, at z s = 2 and A9 = 0.1", the 
values of P(A9) in terms of the Sheth & Tormen (1999) 
mass function are 1.2 times smaller than those derived 
from the PS approach, and this ratio is insensitive to cos- 
mological models. The distribution of the differential lens- 
ing probability normalized by P(A9) based on the Sheth 
& Tormen (1999) mass function remains roughly the same 
as in Figure 1. 

Thirdly, the theoretically predicted distribution of im- 
age separations depends critically on the choice of the cos- 
mological parameters (e.g. Li & Ostriker 2000). Therefore, 
a robust constraint on the WDM particle mass from the 
study of lensing probability for multiple images with small 
angular separations needs a priori precise determinations 
of many other parameters such as CIm, ^a, o~g, and Hq. 

Finally, one may worry about the contamination of the 
supermassive black holes as lenses in the actual applica- 
tion of the gravitational lensing probability to the observa- 
tional tests of WDM model. Indeed, it is believed that all 
the AGNs contain supermassive black holes with masses 
M ~ 10 8 - 10 9 Mq in their centers. The lensing effect 
by these supermassive black holes leads to an image sep- 
aration of typically ~ 0."1 for a quasar at cosmological 
distance, which is indeed of the same order of magnitude 
as the effect produced by small galaxies discussed in the 
present Letter. Yet, the lensing probability of a back- 
ground quasar due to the supermassive black holes resid- 
ing in AGNs is considerably low (~ 10 -8 ) (Augusto & 
Wilkinson 2001). 

Overall, we feel that gravitational lensing should serve 
as a critical test for WDM model although further work 
will be needed to address in more detail various uncertain- 
ties in the issue. 
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Fig. 1. — The differential probabilities of image separations normalized by the total probabilities P(A9) for distant quasars at z B = 1, 2 
and 3 (from bottom to top) in CDM and WDM models. 
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